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ABSTRACT 

We present Spitzer Space Telescope observations of the well-studied extremely red 
objects (EROs) HR 10 and LBDS 53W091 from 3.6 /jm to 160 /mi. These galaxies are 
the prototypes of the two primary classes of EROs: dusty starbursts and old, evolved 
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' galaxies, respectively. Both galaxies, as well as LBDS 53W069, another example of an 
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t~^. ■ HR 10 is detected in the far-infrared. All three EROs have stellar masses of a few 

x 10 11 Mq. Using evolutionary model fits to their multiband photometry, we predict the 
, infrared colors of similar EROs at 1 < z < 2. We find that blueward of observed 10 /mi, 

^ c- 1 the two ERO classes are virtually indistinguishable photometrically. Deep spectroscopy 

and 24 /mi data allow the classes to be separated. 

6 

' Subject headings: infrared: galaxies — galaxies: individual (HR 10) — galaxies: indi- 

ct ■ vidual (LBDS 53W069) — galaxies: individual (LBDS 53W091) 
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1. Introduction 

It has now been nearly twenty years since the first near-infrared surveys identified an intrigu- 
ing population of optically faint galaxies with surprisingly bright near-infrared magnitudes (Elston, 
Rieke, & Rieke 1988, 1989). Initially proposed to be extremely distant [z ~ 10) sources whose op- 
tical emission was absorbed by intergalactic hydrogen (e.g. Mobasher et al. 2005), the "extremely 
red object" (ERO) population 1 is now instead recognized to be primarily comprised of two inter- 
esting galaxy populations: (1) old, evolved galaxies whose red colors are caused by a dearth of 
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1 EROs are typically selected to have extreme optical-minus-infrared colors, such as (R — A")v C ga > 6.0. 
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young, luminous, massive stars, and (2) dusty starburst galaxies whose red colors are caused by 
dust absorption of the bluer light in actively star-forming regions. The prototype of each class 
is, respectively, the "old, dead, and red" radio galaxy LBDS53W091 at z = 1.552 (Dunlop et al. 
1996; Spinrad et al. 1997) and the dusty starburst galaxy HR10 at z = 1.44 (Hu & Ridgway 1994; 
Graham & Dey 1996; Dey et al. 1999). The weak radio galaxy LBDS53W069 at z = 1.432 (Nolan 
et al. 2003; Dey et al. 2006) is a second well-studied example of the old, evolved ERO class. EROs 
have a surface density of approximately 1000 deg -2 for R — K>6 and K < 20 (Vega), and comprise 
~ 10% of all K-selected galaxies to K < 20 (Thompson et al. 1999; Moustakas et al. 2004). For a 
recent comprehensive review of EROs see McCarthy (2004). 

With redshifts 1 < z < 2 (e.g. Cimatti et al. 2003), clustering properties similar to elliptical 
galaxies (Daddi et al. 2003; Moustakas & Somerville 2002), and substantial stellar masses that 
contribute significant fractions of the global stellar mass density at the redshifts they are found 
(e.g. Kong et al. 2006), EROs clearly are an important population in the context of understanding 
galaxy formation and evolution (e.g. Somerville et al. 2004; Nagamine et al. 2005). Volume- 
limited galaxy surveys of the local universe find a clear bimodality in galaxy properties, with the 
"red sequence" generally populated by the more massive and established galaxies (e.g. Kauffmann 
et al. 2003; Blanton et al. 2005). There is increasing evidence that this sequence is securely in 
place not only out to z ~ 1 (Bell et al. 2004; Faber et al. 2006), but possibly to much higher 
redshifts (Rocca-Volmerange et al. 2004; Labbe et al. 2005; Papovich et al. 2006). The processes 
by which these massive galaxies form must happen early. This is consistent with measurements of 
the evolution of the global stellar mass density $l*(z) in field surveys. While between half and three- 
quarters of the present-day stellar mass is in place by z ~ 1 (Dickinson et al. 2003; Rudnick et al. 
2003; Fontana et al. 2004; Drory et al. 2004), only ~ 10% of the present-day stellar mass is in place 
at z > 3 (Drory et al. 2005). The evolution of is dramatic: either the high-redshift accountings 
are incomplete (e.g. because of cosmic variance effects, unaccounted for obscured populations, or 
significant population synthesis deficiencies), or probes in that redshift range will see the most 
dramatic epoch of galaxy assembly in progress. Indeed, studies of infrared selected samples of 
distant galaxies in the Great Observatories Origins Deep Survey (GOODS; Giavalisco et al. 2004) 
find significant evidence of massive galaxy assembly at z ~ 1 — 3 (Caputi et al. 2006; Papovich 
et al. 2006). Understanding ERO demographics and energetics will provide key insight into all of 
these fundamental questions. 

EROs are represented by two different galaxy populations. With enough dedicated time and 
effort, young, dusty starbursts and old, evolved galaxies may be distinguished through spectroscopic 
features in the restframe UV and optical (e.g. Cimatti et al. 2002a). A photometric technique 
proposed by Pozzetti & Mannucci (2000) to distinguish between ERO types using combinations of 
optical and near- infrared colors was designed to work for objects around z ~ 1.5 (e.g. Mannucci 
et al. 2002). Dorman, O'Connell, & Rood (2003) suggest an alternate discriminating technique 
using mid-ultraviolet colors. A more detailed look at the actual A rcst ~ 0.1 — 1.1 /im spectral energy 
distributions (SEDs) of EROs by Moustakas et al. (2004) demonstrates that a larger baseline, 
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extending longward of the observed K-band, is needed to differentiate all but the most dramatic 
EROs. 

In this paper we present 3.6 fim to 160 fim observations of the archetypical EROs obtained with 
the Spitzer Space Telescope (Werner et al. 2004). Section 2 briefly summarizes previous studies of 
these three EROs. Section 3 describes the Spitzer observations and data reduction, followed by an 
analysis of the broadband SEDs of these galaxies in § 4. Surprisingly, the two ERO populations 
remain nearly indistinguishable out to observed 8jum (rest-frame 3.2 ^m) and it is only at longer 
wavelengths that the populations separate (§ 5). We adopt the concordance cosmology and use 
Vega-system magnitudes unless stated otherwise. 

2. Prototypical EROs 
2.1. HR10 

HR 10 (or ERO J164502+4626.4) is the archetype of the dusty starburst subclass of EROs. 
With K' = 18.4 and I — K = 6.0, HR10 was first identified by Hu & Ridgway (1994) in deep, 
multiband imaging of the damped Lya quasar PC 1643+4631A (z = 3.79; Schneider, Schmidt, 
& Gunn 1991). Hu & Ridgway (1994) initially suggested that HR10 and another ERO in the 
field (HR 14 or ERO J164457+4626.0) could be distant, evolved ellipticals lying at z ~ 2 - 3, but 
subsequent observations from the Keck telescopes by Graham & Dey (1996) and Dey et al. (1999) 
showed that HR10 is a moderately distant (z = 1.44) galaxy with an asymmetric morphology and 
[O II] and Ha in emission. This suggested that HR 10 was not related to an evolved elliptical, but 
was rather a distant counterpart of the local ultraluminous infrared galaxies (ULIRGs) discovered by 
the Infrared Astronomical Satellite (IRAS) . Strong detections of HR 10 at submillimeter wavelengths 
(Cimatti et al. 1998; Dey et al. 1999; Greve, Ivison, & Papadopoulos 2003) showed the presence 
of significant dust (Mdust ~ 9 x 10 s M ; Greve et al. 2003), unambiguously demonstrating that 
HR 10 is indeed a gas-rich, dust-enshrouded galaxy. 

A long-standing open issue has been the role of a possible active nucleus in the (long-wavelength) 
energetics of HR 10. Assuming the submillimeter flux is mainly due to optically-thin thermal 
emission from dust heated by a young, star-forming population, the 8 — 1000 /im luminosity, 
Lfir ~ 9 x 10 12 L Q , translates to a star formation rate of about 900 M Q yr _1 (Greve et al. 
2003). This starburst-dominant hypothesis has been favored by the lack of strong evidence of an 
AGN in the spectral data (Dey et al. 1999). Furthermore, the mid-infrared SED of HR10 as 
measured by ISOCAM on the Infrared Space Observatory (ISO) is similar to that of Arp 220, a 
well-studied, local ULIRG whose dominant source of energy is believed to be star formation (Elbaz 
et al. 2002). 

The 850/um photometry of HR 10 remains somewhat uncertain, with several Submillimetre 
Common-User Bolometer Array (SCUBA; Holland et al. 1999) observations resulting in disparate 
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flux density measurements (see Table 1). The observations by Cimatti et al. (1998) and Dey 
et al. (1999) were both obtained using SCUBA in "photometry mode", providing a single pixel 
measurement. The latter data set is twice as large and was obtained in exceptional observing 
conditions and is therefore presumed the more reliable measurement. The observations by Greve 
et al. (2003) were obtained in "jiggle map mode", providing an image of the field and photometry 
closer to that measured by Cimatti et al. (1998). Greve et al. (2003) identify a 6±2 mJy SCUBA 
source approximately 50" SSW of HR 10, and suggest that the extended emission from this second 
source may have affected the sky background estimates in the photometry mode observations, 
thereby corrupting those measurements. An additional source, approximately 30" SE of HR 10, 
may also be confusing the SCUBA measurements. This source is faintly visible in the jiggle map 
of Greve et al. (2003), and is well-detected in our Spitzer observations described in § 3, becoming 
significantly brighter than HR10 at wavelengths longward of 50 /um (see Fig. la). 



2.2. LBDS 53W091 

LBDS53W091 is the archetype of the evolved elliptical subclass of EROs. With K = 18.7 
and R — K = 5.8, LBDS 53W091 is among the reddest counterparts to weak radio sources (1 mJy 
< <Si.4GHz < 50 mJy) in the Leiden-Berkeley Deep Survey (LBDS; Windhorst, van Heerde, & 
Katgert 1984b; Windhorst, Kron, & Koo 1984a). Probing sources with these faint radio fluxes is 
expected to identify sources whose scattered AGN contribution to the total continuum emission 
is small (e.g. Dunlop & Peacock 1993; Eales & Rawlings 1993). Spectroscopic observations with 
the Keck telescope provided an absorption- line redshift, z = 1.552, and a rest- frame ultraviolet 
spectrum very similar to an F6 V star and the local elliptical galaxy M32 (Dunlop et al. 1996; 
Spinrad et al. 1997). The optical spectrum shows no evidence of either an AGN or recent star 
formation. 

The implied stellar age of > 3.5 Gyr was cosmologically-interesting and controversial at the 
time of initial publication (e.g. Spinrad et al. 1997; Krauss 1997; Yi et al. 2000), which was prior 
to the current "concordance cosmology." For example, an S7 m = 1 and Hq = 70kms _1 Mpc -1 
Einstein-de Sitter universe is only 2.3 Gyr old at z = 1.552. Even for the "concordance cosmology" 
(J7m = 1 — f^A = 0.3 and Hq = 70 kms -1 Mpc -1 ), the universe is only 4.1 Gyr old at z = 1.552, 
requiring an early formation of the stellar populations of LBDS 53W091. Galaxy 3a, approximately 
3" SE of LBDS 53W091 and a few tenths of a magnitude fainter, has an optical spectrum and optical 
to near-IR colors similar to LBDS 53W091. Galaxy 3a is thought to be a faint companion, of similar 
age and redshift to LBDS53W091 (Spinrad et al. 1997; Bunker et al. 2006). 

Images of LBDS 53W091 with the Near-Infrared Camera and Multiobject Spectrograph (NIC- 
MOS) on the Hubble Space Telescope (HST) reveal LBDS 53W091 to have a regular r 1 / 4 de Vau- 
couleurs (1948) profile with an effective radius of 0"32±0 / /08 (2.7±0.7 kpc; Waddington et al. 2002; 
Bunker et al. 2006). Assuming passive evolution of its stellar populations, LBDS53W091 lies well 
within the scatter of the Kormendy (1977) relation, consistent with this milliJansky radio source 
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at z ~ 1.5 being an early-type galaxy which will passively evolve into an ordinary L* elliptical by 
the present day. 

2.3. LBDS 53W069 

LBDS 53W069 is a second, well-studied example of the evolved elliptical subclass of EROs. 
Slightly brighter than LBDS 53W091 at optical and near-infrared wavelengths, LBDS 53W069 is 
a weak radio source whose optical spectrum reveals an ps 4.0 Gyr-old stellar population at z = 
1.432 (Nolan et al. 2003; Dey et al. 2006). NICMOS images again reveal a regular, r 1 / 4 profile, 
consistent with the Kormendy relation (Waddington et al. 2002; Bunker et al. 2006). Submillimeter 
observations of LBDS 53W069 obtained by Archibald et al. (2001) fail to detect this galaxy, again 
consistent with LBDS 53W069 no longer being actively involved in star formation. 

3. Observations and Data Reduction 

We obtained images from 3.6 /im to 160 /im of these prototypical EROs with the Spitzer Space 
Telescope as part of guaranteed time observations. These observations employed both the Infrared 
Array Camera (IRAC; Fazio et al. 2004) and the Multiband Imaging Photometer for Spitzer (MIPS; 
Rieke et al. 2004). 

3.1. IRAC 

Both LBDS 53W091 and LBDS 53W069 were observed with IRAC on 2004 February 17, and 
HR 10 was observed on 2004 February 18. Exposures with IRAC totaled 900 sec per target, obtained 
as nine 100 sec dithered frames using the small-scale, cycling dither pattern. Due to the higher 
background at longer wavelengths, there are twice as many frames in IRAC channel 4 (8.0 /im), 
each with half the frame time. 

We made minor modifications to the pipeline-processed data. Beginning with the basic cali- 
brated data from the IRAC pipeline (version Sll.0.2), we corrected the "pulldown" and "muxbleed" 
associated with well-exposed pixels in IRAC channels 1 (3.6 /j,m) and 2 (4.5 //m). The empirically- 
derived muxbleed corrections were derived from a large set of IRAC data taken with the same 
exposure time. We then reran the MOPEX 2 mosaicing software with modified cosmic-ray rejection 
parameters and an output pixel scale reduced by 50% to improve the sampling of the final images. 
Source Extractor (Bertin & Arnouts 1996) was used to derive photometry in 5" diameter apertures 
for most sources, which were then corrected to total magnitudes assuming sources are unresolved 
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at the relatively coarse (~ 2") spatial resolution of Spitzer. The exceptions were galaxies in close 
proximity to other sources, e.g. LBDS53W091. For these galaxies, where source blending can 
compromise accurate photometry, we measured fluxes in 2" diameter apertures, which were then 
corrected to total magnitudes assuming the sources were unresolved. 

Fig. 1 presents the processed IRAC images. For reference, optical images obtained with the 
Wide-Field Planetary Camera 2 (WFPC2; Trauger, Ballester, & Burrows 1994) on the HST are 
also presented. These data, all obtained with the F814W filter, are from the Multimission Archive 
at the Space Telescope Science Institute 3 . 

3.2. MIPS 

Both HR10 and LBDS 53W069 were observed with MIPS on 2004 April 8, and LBDS 53W091 
was observed on 2004 April 10. The observations were performed in MIPS small scale photometry 
mode with 30 sec x 5 cycles, 10 sec x 10 cycles and 10 sec x 5 cycles at 24, 70 and 160 /im, 
respectively. The 160 /tin observations also incorporated a 1 x 3 map, stepping by 75% of the field 
of view. The total integration times on source were ~ 2100 sec at 24 /xm, 1000 sec at 70 //m and 
300 sec at 160/im. 

The 24 /im and 70 /im pipeline (version SI 1.4.0) basic calibrated data (BCD) products were the 
starting point for the reductions. The 24 /im BCDs were then reprocessed to correct for flat fielding 
patterns due to residue on the MIPS scan mirror, "jailbar" patterns caused by bright cosmic rays 
artificially suppressing one or more of the four readout channels, and gradients in the background 
subtraction due to variation in either the dark current or the actual sky background. The 70 /im 
BCDs were column and time filtered to correct for cosmic rays and response variations. Details of 
24 /im and 70 //m artifacts can be found in the MIPS data handbook. The reprocessed BCD frames 
were mosaiced together using MOPEX. For 160 /im, the default pipeline processed mosaic failed to 
detect any of the sources. For this bandpass, we present 3cr upper limits derived from the Spitzer 
Sensitivity Performance Estimation Tool 4 for medium background level. This is a conservative limit 
since all three sources were actually observed at low background levels (S^m < ^ MJy sr _1 ). 

Photometry at 24 /im and 70 /im was performed using a point source fitting algorithm, whereby 
the positions of the source in the shorter wavelength IRAC data was fed in as a prior. This has 
been shown to provide accurate photometry as demonstrated in the GOODS datasets (Chary et 
al., in preparation). The 70 /im photometry derived here was corrected upwards by approximately 
5% to account for flux suppression by the time and column filtering process. 

The Spitzer photometry is presented in Tables 1 — 3, along with photometry at other wave- 
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lengths derived from the literature. Where no detection is made, 3a upper limits are given. The 
broadband SEDs of all three EROs are presented in Fig. 2. 

4. Spectral Energy Distributions 

We next fit the broadband SEDs with galaxy template models with the goal of deriving basic 
properties of the galaxies, such as stellar mass and far- infrared luminosity. The new data from 
Spitzer have the capacity to significantly refine previous determinations of these properties. For 
the SED fitting, described below, we omit some of the photometric points presented in Table 1. 
First, where HST imaging is available at similar wavelengths to ground-based imaging (e.g. i- 
band/F814W and J-band/FHOW), we use only the higher signal-to-noise ratio HST photometry. 
Second, we do not include < 3a detections in our fits. Finally, we omit the early (circa 1985) CCD 
Palomar 200" gri imaging of the LBDS sources whose photometry typically compare poorly with 
more recent, deeper photometry obtained by HST and the Keck telescopes. 

We fit Bruzual & Chariot (2003, BC03) model SEDs to the optical thru mid-infrared (A D bs < 
8/xm) fluxes of these sources. We adopt a Salpeter initial mass function with solar metallicity 
templates, leaving the mass of the galaxy, dust extinction, age of the starburst and e-folding 
time of the starburst as free parameters. Dust was assumed to be distributed in a screen. Both 
starburst and galactic extinction laws were tried. The models were redshifted and convolved through 
the relevant filter response/quantum efficiency curves. A weighted, reduced, chi-square fit was 
performed to the observed data points after multiplication with a model normalization factor which 
represents the combined mass of stars and gas. The mass fraction in stars was extracted from the 
corresponding BC03 model file and multiplied by the normalization factor to derive the total mass 
in stars. Fig. 2 presents the best-fit BC03 models to all three EROs and Table 4 presents the 
parameters of the best-fit models. Fig. 2a presents the full SED, while Fig. 2b only shows the bluer 
wavelengths of the SED, where the emission is dominated by starlight. 

We find that the LBDS galaxies have stellar masses of ~ 3 x 10 11 M Q , while HR 10 has a stellar 
mass roughly twice as large. Modulo differences in the assumed initial mass function, these values 
are 4 — 8 x the average stellar mass of an L* galaxy in the local universe (M* ~ 7 x 1O 1O M0; 
Cole et al. 2001) or of evolved 2 < z < 3 infrared-selected galaxies in the Hubble Deep Field 
- South (M* ~ 8 x 1O 1O M ; Labbe et al. 2005). We find ages of a few Gyr for the evolved 
galaxies, comparable to the more carefully determined ages based on fitting the observed optical 
spectra (Spinrad et al. 1997; Nolan et al. 2003). Surprisingly, our best-fit model to HR10 has a 
comparable age, albeit with heavy extinction, Ay = 2.4. 

To fit the longer wavelength Spitzer and (sub)mm data for HR 10, we adopted the Chary h 
Elbaz (2001, CE01) range of far-infrared spectral energy distributions. The CE01 templates pro- 
vide a range of mid- and far-infrared SEDs as a function of the far-infrared luminosity, Lfir = 
L(8 — 1000/um). They are empirically derived templates based on observations of nearby infrared 
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luminous galaxies and incorporate components associated with continua from multi-temperature 
dust grains, emission from poly cyclic aromatic hydrocarbons, and absorption due to silicate dust. 
These templates have been shown to follow the luminosity-luminosity correlations in IRAS pass- 
bands seen for nearby galaxies. As with the BC03 fits to the stellar emission, the dust emission has 
been redshifted and convolved through the relevant passbands. We also attempted to fit the ob- 
served fluxes with a combination of the CE01 templates and obscuration based on the Li & Draine 
(2002) model which would most prominently introduce a spectral feature at 9.7/im corresponding 
to silicate absorption. We find that a fit including the Li & Draine (2002) model for dust obscura- 
tion provides a better fit to the mid-infrared data points, although the derived -Lfir = 3 x 1O 12 L 
remains essentially unchanged. The resultant obscuration corresponds to Ay ~ 11 mag, which is 
significantly larger than that seen in the optical/near-infrared. This suggests that HR10 has regions 
of heavily obscured star- formation which make it a ULIRG. With more accurate measurements of 
the Wien side of the thermal dust emission, our derived far-infrared luminosity is approximately 
40% that derived by Greve et al. (2003) and Elbaz et al. (2002). Fig. 2a presents the best-fit CE01 
far-infrared template. Since the CE01 templates are based on real galaxies and include stellar 
populations, Fig. 2a shows the best-fit CE01 template only at longer wavelengths (A rcst > 5 |im), 
stitched to the best-fit BC03 stellar population template at shorter wavelengths. 

It has been suggested that HR 10 is a scaled up version of Arp 220 (Elbaz et al. 2002). We find 
that the mid-infrared and far- infrared SED of HR 10 shows significant differences in comparison 
with Arp 220 (Fig. 3). The mid-infrared colors of HR 10 are better fit by an SED including strong 
polycyclic aromatic hydrocarbon (PAH) in emission and weak silicate in absorption. The derived 
extinction from the silicate absorption is Ay ~ 10 mag. In contrast, Arp 220 has a very strong 
silicate feature, indicative of Ay > 50 — 80 mags. Since HR 10 is approximately twice as luminous 
than Arp 220 (Elbaz et al. 2002), the strength of the silicate feature apparently does not scale with 
far-infrared luminosity as would be naively expected. HR10 also appears to have a cooler color 
temperature for the far-infrared emission which is more typical for the local ULIRG population 
than the SED exhibited by Arp 220. 

Since the photometry of HR 10 is well-fit with purely stellar populations with dust, we infer 
that the energetics of this galaxy is dominated by starlight. The radio detection of HR 10 is fully 
consistent with the CE01 template, implying that it obeys the radio - far-infrared correlation 
observed in local star- forming galaxies (Yun & Carilli 2002). If HR10 does indeed host an AGN, 
the active nucleus apparently plays a negligible role in the broadband emission from this galaxy. 

Finally, to fit the radio portion of the LBDS SEDs, we fit the data with a self-absorbed 
synchrotron emission spectrum. Neglecting the high-energy cut-off, this can be approximated by 
a parabola-shaped curve in the log-log plane for frequencies lower than the frequency at which the 
plasma optical depth becomes unity, v t (e.g. Andreani et al. 2002). We adopt v t = 31.6 MHz for 
both sources and fit for the normalization and curvature of the radio spectrum. The results are 
plotted in Fig. 2a and illustrate that the synchrotron emission which dominates the radio portion 
of the LBDS SEDs has negligible contribution to the optical through far-infrared emission for these 
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sources. 

5. Separating Dusty Starbursts from Old Ellipticals 

It has been nearly a decade since the recognition that EROs constitute a bimodal population of 
dusty, star-forming galaxies and quiescent, evolved galaxies. As EROs seem to be the likely z ~ 1.5 
progenitors of present day early- type galaxies (e.g. Moustakas & Somerville 2002), measuring their 
relative fractions with confidence is an important measurement for feedback into models of galaxy 
formation. 

Deep spectroscopy can unambiguously separate the two ERO populations, but is very expensive 
in telescope resources and typically suffers from high incompleteness rates. Independent spectro- 
scopic samples of « 25 EROs studied by Cimatti et al. (2002b) and Yan, Thompson, & Soifer 
(2004a) both found roughly equal distributions of emission line galaxies and purely absorption line 
galaxies. These spectroscopic classes are presumed associated with dusty starbursts and evolved 
early-type EROs, respectively. More recently, Doherty et al. (2005) reported on a spectroscopic 
sample of EROs nearly twice as large as the previous studies, finding that 75% of their sample have 
a dominant old stellar population, but only 28% of the sample is completely devoid of evidence of 
recent star formation. With typically one-third of the observed EROs lacking sufficient signal to 
enter the spectroscopic samples, the error bars on these population statistics are substantial. 

In principle, high resolution morphologies from EST would be useful at distinguishing dusty 
starbursts from quiescent, old galaxies. This expectation is supported by Fig. 1, which shows that 
the dusty starburst HR 10 has a disturbed morphology at observed /-band, while both LBDS 53W091 
and LBDS 53W069 have regular, undisturbed morphologies. From a deep HST WFPC2 archival 
study of 115 EROs, Yan & Thompson (2003) found that 30 ± 5% of EROs are bulge-dominated, 
64 ± 7% are disk-dominated, and only 6% defy morphological classification. However, from their 
spectroscopic study, Yan et al. (2004a) found very little correlation between spectroscopic classifica- 
tion and morphological classification, with both emission line galaxies and absorption line galaxies 
evenly comprised of bulge- and disk-dominated galaxies. Similarly, from a large, magnitude-limited 
sample of EROs in the GOODS-South field, Moustakas et al. (2004) found that the mean SEDs 
at 0.3 — 2.2 / um are basically independent of morphology. The implication is that there is no robust 
way to distinguish ERO subclasses based purely on imaging data shortward of ps 3^m. 

However, Spitzer mid-infrared data enable the ERO subclasses to be separated. Fig. 4 presents 
the best-fit ERO SEDs (§ 4) redshifted to 1 < z < 2, the range at which most EROs are thought 
to reside. Despite previous claims that J — K might be a useful discriminant between the ERO 
subclasses (e.g. Pozzetti & Mannucci 2000), Fig. 4 clearly demonstrates that ground-based near- 
infrared data cannot robustly distinguish the subclasses. 

In theory, the best-fit SEDs suggest that IRAC data out to 8/xm might be useful at separating 
the populations. While the quiescent, evolved galaxies have SEDs which are falling longward of 
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the 1.6/xm stellar peak, the star-forming SED rises into the mid- and far- infrared due to the cold 
gas and dust associated with star formation. Though the models differ by 1 — 1.5 mag in K— [8.0] 
and [3.6] — [8.0] colors, actual observations show more subtle color differences. Not surprisingly, 
the 100 x difference in the 24 /Ltm properties of the best-fit SED provide the most robust separation 
of the star-forming and evolved galaxy EROs. 

6. Conclusions 

We present Spitzer IRAC and MIPS observations of three well-studied EROs. HR10 is a 
dusty starburst class of EROs, while the weak radio sources LBDS 53W091 and LBDS 53W069 
are evolved, quiescent EROs. The IRAC data show the 1.6/um stellar peak in all three sources, 
providing accurate derivations of basic properties of their stellar populations. Consistent with 
the bright near-infrared magnitudes, we find all three galaxies are massive, with stellar masses 
M* = (3 - 6) x lO n M . We also derive ages of a few Gyr for all three galaxies, though HR 10 is 
found to be dustier and to have proportionally more recent star formation. 

Only HR 10 is detected by MIPS, illustrating that 24^m photometry provides a robust dis- 
criminant between the two primary classes of EROs. Blueward of 10/xm, we find that both ERO 
classes are similar; identifying the relative fractions of massive galaxies at z ~ 1.5 which are old or 
dusty and starforming requires either longer wavelength data or deep spectroscopy. From a sample 
of EROs in a 64 square arcmin area of the Spitzer First Look Survey, Yan et al. (2004b) find that 
approximately half of the ERO population is detected to a 24/im flux limit of 40/uJy (3a), suggest- 
ing a roughly even distribution of starbursts and quiescent galaxies in the ERO classification. This 
breakdown is consistent with the courageous, yet incomplete, spectroscopic studies (e.g. Cimatti 
et al. 2002a; Yan et al. 2004a; Doherty et al. 2005). Since EROs are the likely progenitors of 
massive early-type galaxies in the early universe, this suggests both that substantial formation of 
early-type galaxies occurs at z > 2, but that many massive galaxies are also still actively forming 
stars at z < 2. 
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Table 1. Photometry of HR 10 



r\i .. . ] 

Observed 


Rest 




Detector/ 




Wavelength 


Wavelength 


Flux Density 


Instrument 


Reference 


3620 A(J7) 


1480 A 


< 0.03 /iJy 


WHT 


1 


4400 k(B) 


1800 A 


0.16 ± 0.07 /iJy 


UH 2.2 m 


2,3 


4900 A(G) 


2010 A 


0.09 ±0.04 /iJy 


WHT 


1 


5420 A(V) 


2220 A 


< 14 » Tv 


WHT 


1 


6450 A(R) 


2640 A 


21 ± 08 /yTv 


WHT 


1 


8140 A(F814W) t 


3340 A 


0.52 ± 0.06 /uJy 


WFPC2/HST 


4 


8200 A(J) 


3360 A 


0.58 ± 0.22 /xJy 


WHT 


1 


8490 A(7) 


3480 A 


0.41 ± 0.13 /uJy 


UH 2.2 m 


2 3 


1.2 fim (J) f 


4920 A 


6.4 ± 2.1 /iJy 


UH 2.2 m 


2 3 


1.6 fim (H)i 


6560 A 


14.8 ± 3.6 /iJy 


UH 2.2 m 


2 3 


2.2 (Um (_K")t 


9020 A 


27.7 ± 0.6 /uJy 


Keck/NIRC 


2 


3.6 ftm 1 


1.5 fim 


65.06 ± 6.53 /iJy 


IRAC / Spitzer 


5 


4.5 ftm^ 


1.8 //in 


80.76 ± 8.11 /iJy 


IRAC / Spitzer 


5 


5.8 ftm 1 


2.4 |im 


59.24 ± 6.41 /xJy 


IRAC / Spitzer 


5 


8.0 /xmt 


3.3 fim 


57.43 ± 5.96 /xJy 


IRAC / Spitzer 


5 


12 /zmt 


4.9 /xm 


85 ± 50 /uJy 


ISOCAM/ZSO* 


6 


15 /im^ 


6.1 /xm 


203 ± 62 /iJy 


ISOCAM/ZSO* 


6 


24 /mi ' 


9 8 nm 


350 + 50 /i.Tv 


MIPS 1 Snitzer 


5 


70 //m f 


28.7 /im 


5.6 ± 2.0 mJy 


MIPS /Spitzer 


5 


90 /im 


36.9 /im 


< 120 mJy 


ISOPHOT/ISO 


6 


100 jum 


41 /im 


< 270 mJy§ 


IRAS 


4 


160 /^m* 


66 /xm 


< 50 mJy 


MIPS/ Spitzer 


5 


170 /im 


70 /im 


< 120 mJy 


ISOPHOT/JS'O 


6 


450 /um'f 


184 /im 


32.2 ±8.5 mJy 


SCUBA/ JCMT 


4 


850 fim 


348 /im 


8.7 ± 1.6 mJy 


SCUBA/JCMT 


7 


850 /im 


348 /im 


4.89 ± 0.74 mJy 


SCUBA/JCMT 


4 


850 yum 1 


348 /im 


8 ± 2 mJy 


SCUBA/JCMT 


8 


1350 /an* 


553 /im 


2.13 ±0.63 mJy 


SCUBA/JCMT 


4 


3.6 cm* 


1.5 cm 


35 ± 11 /iJy 




9 


21 cm 


8.2 cm 


< 300 /iJy § 




4 



^Data used in SED modeling (§4). 

■'■Systematic and statistical errors have been added in quadrature. 
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§ Reference does not specify significance of non-detection. 

References. — (1) Haynes et al. (2002), (2) Graham & Dcy (1996), (3) Hu & Ridgway (1994), (4) Dey 
ct al. (1999), (5) This paper, (6) Elbaz et al. (2002), (7) Cimatti et al. (1998), (8) Greve et al. (2003), (9) 
Frayer (1996) 
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Table 2. Photometry of LBDS 53W091 



Observed 


r? Q pf 




Detector/ 




Wavelength 


Wavelength 


Flux Density 


Instrument 


Reference 


4960 A(g) 


1940 A 


0.16 ± 0.07 jiJy 


Palomar 200" 


1 


6420 A(i?)t 


2520 A 


0.47 ± 0.09 jiJy 


Keck/LRIS 


2 


6480 A(r) 


2540 A 


0.14 ±0.06 jiJy 


Palomar 200" 


1 


8140 A(F814W) t 


3190 A 


0.80 ± 0.02 jiJy 


WFPC2/iZST 


3 


8160 A(i) 


3200 A 


0.45 ± 0.08 jiJy 


Palomar 200" 


1 


1.1 (F110W) t 


4310 A 


3.44 ± 0.10 jiJy 


NICMOS/ HST 


3 


1.2 /jm (J) 


4700 A 


10.00 ± 0.92 j/Jy 


UKIRT 


1 


1.6 /xm (fT)t 


6270 A 


15.85 ± 1.46 j/Jy 


UKIRT 


1 


2.2 jan (if)t 


8620 A 


21.48 ± 2.58 j/Jy 


UKIRT 


1 


3.6 /imt 


1.4 


35.63 ± 3.64 /xJy 


IRAC /Spitzer 


4 


4.5 fim' 


1.8 fim 


37.22 ± 3.82 /iJy 


IRAC/ Spitzer 


4 


5.8 /imt 


2.3 /xm 


19.98 ± 3.01 /iJy 


ITLAC/ Spitzer 


4 


8.0 /imt 


3.1 /7,rn 


14.23 ± 2.04 /iJy 


IRAC /Spitzer 


4 


24 /vm' 


9.4 /7,m 


< 40 /xJy 


MIPS /Spitzer 


4 


70 /viii' 


27 /7,m 


< 8.1 mJy 


MIPS/ Spitzer 


4 


160 /jm* 


63 /xm 


< 50 mJy 


MIPS/ Spitzer 


4 


6.2 cm^ 


2.4 cm 


6.5 ± 0.4 mJy 




2 


19 cmt 


7.4 cm 


23.0 ± 1.7 mJy 




2 


21 cmt 


8.2 cm 


22.1 ±2.0 mJy 




1 


50 cm' 


20 cm 


66.0 ± 3.9 mJy 




1 



"'"Data used in SED modeling (§4). 



References. — (1) Waddington et al. (2000), (2) Spinrad et al. (1997), (3) Waddington et al. 
(2002), (4) This paper 
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Table 3. Photometry of LBDS 53W069 



Observed 


Rest 




Detector/ 




Wavelength 


Wavelength 


Flux Density 


Instrument 


Reference 


4960 k(g) 


2040 A 


0.10 ±0.06 /xJy 


Palomar 200" 


1 


6480 A(r) 


2660 A 


0.30 ± 0.06 /xJy 


Palomar 200" 


1 


8140 A(F814W) t 


3350 A 


0.77 ±0.02 /xJy 


WFPC2/#ST 


2 


8160 A(i) 


3360 A 


0.46 ± 0.06 /xJy 


Palomar 200" 


1 


1.1 /xm (FllOW)t 


4520 A 


3.80 ± 0.07 /xJy 


NICMOS/ HST 


2 


i.z /im (J ) 


A OQfl A 

4yoU A 


iz.oy zt i.oo /iJy 


Urvlrtl 


i 
i 


1.6 /xm (fl-)t 


6580 A 


12.59 ±2.10 /xJy 


UKIRT 


1 


2.2 /xm (K)t 


9050 A 


24.43 ± 2.48 /xJy 


UKIRT 


1 


3.6 /imt 


1.5 fim 


43.48 ± 4.38 /xJy 


IRAC /Spitzer 


3 


4.5 /im^ 


1.9 fim 


43.13 ± 4.35 /iJy 


IRAC /Spitzer 


3 


5.8 /irat 


2.4 /7,rn 


20.72 ± 2.92 fiJy 


IRAC /Spitzer 


3 


8.0 /imt 


3.3 /im 


22.18 ± 2.70 fiJy 


IRAC / Spitzer 


3 


24 /7,m' 


9.9 /im 


< 50 /iJy 


MIPS/ Spitzer 


3 


70 /xm' 


29 /im 


< 7.8 mJy 


MIPS/ Spitzer 


3 


160 /xmt 


66 /im 


< 50 mJy 


MIPS/ Spitzer 


3 


450 /xm 


185 /jm 


< 49 mJy 


SCUBA/JCMT 


4 


850 /xm 


360 /xm 


< 3.12 mJy 


SCUBA/JCMT 


4 


21 rant 


8.6 cm 


3.7 ±0.3 mJy 




1 


50 rant 


21 cm 


7.8 ± 0.9 mJy 




1 



tData used in SED modeling (§4). 



References. — (1) Waddington et al. (2000), (2) Waddington et al. (2002), (3) This paper, (4) 
Archibald et al. (2001) 
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Table 4. BC03 Model Fits to ERO Stellar Populations 





HR10 


LBDS 53W091 


LBDS 53W069 


M* [10 n M Q ] 


5.6 


3.0 


3.2 


Ay [mag] 


2.4 


0.4 


1.4 


age [Gyr] 


3.0 


3.5 


4.0 


e-folding [Gyr] 


2.0 


0.7 


2.0 


^FIR [Lq] 


3 x 10 12 







HST/WFPC2 - F814W Spitzer/IRAC - [5.6] Spitzer/IRAC - [4.5] 




20 10 -10 -20 20 10 -10 -20 20 10 -10 -20 



Spitzer/IRAC - [8.0] Spftzer/MIPS - [24] Spitzer/MIPS - [70] 




20 10 -10 -20 20 10 -10 -20 20 10 -10 -20 

Arcsec Arcsec Arcsec 



Fig. la. — HST and Spitzer images of the prototypical ERO HR10, located at a = 
16 h 45 m 02.36 s ,5 = +46°26 / 25 / /5 (J2000). Cameras and filters are indicated. The fields of view 
shown are 50" on a side and relevant sources are labeled. North is up, and east is to the left. The 
inset is 4" on side. Note the galaxy approximately 30" SE of HR 10 which dominates at 70 /urn. This 
source also dominates the 160 ^irn image; we suggest that it may be contaminating the published 
850 /urn photometry as well. 



-19- 



HST/WFPC2 - F814W Spitzer/IRAC - [5.6] Spitzer/IRAC - [4.5] 




20 10 -10 -20 20 10 -10 -20 20 10 -10 -20 




Fig. lb. — HST and Spitzer images of the prototypical ERO LBDS53W091, located at a = 
17 h 22 m 32.63 s ,(5 = +50°06'0l'. / 5 (J2000). Cameras and filters are indicated. The fields of view 
shown are 50" on a side and relevant sources are labeled. North is up, and east is to the left. The 
inset is 4" on side. LBDS 53W091 is not detected in the MIPS bands. 
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Fig. lc. — HST and Spitzer images of the prototypical ERO LBDS 53W069, located at a = 
17 h 20 m 02.5 s ,5 = +49°44'51" (J2000). Cameras and filters are indicated. The fields of view shown 
are 50" on a side. North is up, and east is to the left. The inset is 4" on side. LBDS 53W069 is not 
detected in the MIPS bands. 
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Restframe Wavelength for z=1.5 (/xm) 
0.1 1 10 1 00 1 000 1 4 1 5 

= i i iiij 1 — i i 1 1 1 n| 1 — i i 1 1 1 1 1| 1 — i i 1 1 1 1 1| 1 — i i 1 1 1 iij 1 — i i 1 1 1 1 1 1 1 — i i 1 1 1 1 1| 1 — n 




0.1 1 10 100 1000 10 4 10 5 10 6 



Observed Wavelength (yu-m) 

Fig. 2a.— SEDs of HR10, LBDS53W091, and LBDS 53W069. The topmost axis indicates rest- 
frame wavelengths for z = 1.5, the approximate redshift of all three galaxies. Photometry not used 
in the SED fits are indicated with open symbols (see §4). The solid line shows the best-fit SEDs. 
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Restframe Wavelength for z=1.5 (/xm) 



0.1 1 




1 10 



Observed Wavelength (/J,m) 

Fig. 2b.— SEDs of HR10, LBDS53W091, and LBDS 53W069, restricted to the stellar portion 
of the spectrum. The topmost axis indicates restframe wavelengths for z = 1.5, the approximate 
redshift of all three galaxies. Photometry not used in the SED fits are indicated with open symbols 
(see §4). The solid line shows the best-fit SEDs. 
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Fig. 3. — Mid-infrared spectrum of Arp 220, shifted to z = 1.44 and scaled to approximate the 
A < 20/um photometry for HR10 (filled circles). HR10 has significantly less silicate absorption 
than Arp 220 and appears to have cooler dust dominating the far-infrared emission. 
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• HR10 

■ LBDS 53W091 

□ LBDS 53W069 



_L 



_L 



_L 



1.2 



1.4 1.6 
Redshift 



1.8 




1.4 1.6 
Redshift 




Fig. 4. — Infrared colors (in the AB system) for the best-fit SEDs redshifted to 1 < z < 2, the 
expected redshift range of EROs. A Rayleigh-Jeans tail has been assumed for the LBDS galaxies 
at A bs > 25/im. Colors of the prototypical EROs are also shown with uniform 20% error bars. 
Though the models suggest that IRAC 8/xm data might be effective at separating dusty starbursts 
(e.g., HR10) from evolved stellar populations (e.g., LBDS53W091 and LBDS53W069), actual 
photometry implies that MIPS 24/im data is necessary. 



